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Abstract: We present a new parametric class of spherically symmetric analytic solutions of 
the general relativistic field equations in canonical coordinates, which corresponds to causal 
models of perfect fluid balls. These solutions describe perfect fluid balls with infinite central 
pressure and infinite central density though their ratio is positively finite and less then one. 
From the solutions of this class we have constructed two causal models in which outmarch of 
pressure, density is positive and monotonically decreasing and pressure -density ratio is less 
than one throughout with in the balls. Corresponding to these models we have maximized the 
Neutron star masses 3.24M and 3.48M @ with the linear dimensions 32.09Kms and 34.36Kms 
respectively with equal surface red shift 0.5811. 

1 .Introduction 

Numerous attempts have been reported to obtain parametric classes of exact 
solutions of Einstein's field equations representing perfect fluid ball in 
equilibrium. [1] -[17] These solutions have four arbitrary constants. The usual 
boundary conditions determine three arbitrary constants leaving one undetermined 
such a solution represents a class of solutions, undetermined constant being a 
parameter. The undetermined constant to be so adjusted to obtain suitable equation of 
state i. e. positive and monotonically decreasing expressions for pressure and density, 
pressure-density ratio must be less then 1 and principle of causality should be obeyed 
within the ball. In this paper we present yet another parametric class of solutions. 

2.Field Equations and Method of Obtaining Analytic Solutions. 

We consider the static and spherically symmetric metric in canonical co-ordinates 

ds 2 = -e x dr 2 -r\d6 2 +sin 2 6d<f> 2 ) + e v dt 2 (1) 

where A, andf are functions of r .Einstein's field equations of gravitation for a non 

empty space- time are 
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where R tj is a Ricci tensor , TV is energy- momentum tensor and R the scalar curvature. 
The energy-momentum tensor TV is defined as 

7;. = {p+pc 2 )v i v J -pg v (3) 
where p denotes the pressure distribution, p the density distribution and v j 
the velocity vector, satisfying the relation 

g,vV=l (4) 
Since the field is static, therefore, 

v 1 = v 2 = v 3 = and v 4 = ^=L= (5) 

Thus we find that for the metric (1) under these conditions the field equation 
(2) reduces the following: 

871G _ x u' 11 

— ^ = e (— + —)-— (6) 
c r r r 
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— P = e (~ — + — + — ) (V) 
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— p = e 2 ( r ) + — (8) 

c r r r 

where prime (/) denotes differentiation with respect to r . From equations (6) and (7) 
Tolman [1] obtained following differential equation in A and u . 



<ir 2r 



2r 



(9) 



By assigning one of the field variables A and t> as some known functions of r the 
equation (9) reduces into a form which on integration defines the metric(l) 
completely and the fluid parameters p and p an be calculated from equation (6) and 
equation (8). Many attempts have been made to solve equation (9) by assuming an 
adhoc relationship between v and A, [2], [4], [9] . 

In this paper we use a substitution. 



e /2 = U, V = e 



(10) 



Thus equation (9) reduces to the following Linear differential equation in V. 



dV 
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dr dr 
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On solving (1 1) we get, 



(11) 



V = e~ l = 



(rU'+uf 



. .Aril 1 +u)uJ^^dr 
A ~ 1 \ p 



AUdr 



r AUdr 
J r{rU'+u) 



(12) 



Where A is an arbitrary constant. Our task is to explore the possibilities of choosing 
U such that the right hand side of equation (12) becomes integrable. In this paper we 
assume 

AUdr 

(13) 



r AUdr 



I and n are arbitrary constants. Equation (13) results into a second degree 
homogenous differential equation in U. 



iirL + 2n)rU' +(l-4p = 



(14) 



The solution is 



where 
n — l 

a 



2n 2n 



U = C,r a+b - 1 +C 2 r a - h - 1 (15) 



b = —4{n-lf+\6n, (16) 



(17) 



Provided, n ^ . 
C x and C 2 are arbitrary constants. Also(12) is simplified into 

v _-x_ r &+ "-'- (a - h){ " +2) [A-2l] 

~\a + b)cy h +{a-b)C 1 p 
where 

/ = \r l - n - 9+ ^ n+2 \a + b)C x r 2b + {a -b]C 2 '\ + '[c i r 2b + C 2 \lr (18) 

The solution is complete if (18) is integrated. In the foregoing sections we shall study 
a method of solving (18) and present a detail discussion of the resulting solution. It 
may be mentioned here that for n =0,we rediscover the class of solutions given by 
Tolman (Tolman's V solutions). 



3. New Class of Solutions 

The equation(18) can be integrated by the method of substitution, if we assume 

2b-\ = l-n-9 + {a-b\n + l) (19) 
In view of (16), the equation (19) yields a quadratic equation in / 



3(n + l)/ 2 +2(n 2 -Un)l + 6n 3 -\3n 2 +57n = (19a) 
which solves into 



-n{n-\\) + J-n{\ln +n 2 +1 In + 171) 

/ = — ^ V ~ , if (n+l)+ ( 20) 

3(n + 1) 

/ = 3.1667 ,i{(n+l)=0 (20a) 



Here, we have considered only the positive radical sign ,as corresponding to negative 
radical sign , e v becomes singular at the origin. We thus obtain a new class of 
solutions of equation (10) as follows: 



f = ( C ir a+h - 1 +C 2 r a - h - 1 f (21) 



Ar~ 2b 



e _x = T t-^-b(cJ +C'r- 2b ) (22) 

[{a + b)cy b +(a-b)C 2 Y 2 1 ' 

where 

C/ = (a + + 2 )C, , C 2 ' = [ a(n + 2)+(n + 4)b ]C 2 (23) 



5 = "7 xu v r~ ( 24 > 

b(a + bf(n + 2\n + 3)C, 

We observe that ( e u ) r=0 becomes singular for all values except - 1 .3684 < n < -1 and 
- 0.03 1 < n < . It may be pointed out here that the class of solutions obtained by 
Wyman M (Wyman 1949) has expression for e"is similar to (21), however, the two 
classes of solutions are disjoint. 



4. Properties of the New Class of Solutions 

In view of (21) and (22) we obtain from (6) and (8), the pressure and density 
distribution respectively. 



Ar 



8xG 1 



[(a + bfcr 2 " +{a-b)C 2 f 
' (2a + 2b - l)C, + (2a -2b-\)C 2 r 



*(c/+C 2 'r-») 
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C x +C 2 r 



-2b 



-1 



(25) 



— P = 



A[(a + b){(n + 3)26 - \}C X + (2b - l)(a - Z))C 2 r' 2fc ] | 
' [(a + b)C x r 2h +(a-b)C 2 ] n+ ' + 

+ (l-2b)BC 2 r- 2h +1 + 5C/ 



(26) 



In addition to the parameter n ,the solutions (21) and(22)contain three arbitrary 
constants Cj, C2 and A. These are to be determined by matching the solutions (21) 
and (22) with Schwarzschild exterior solution for a ball of mass M and linear 
dimension 2r h : 



p(r b ) = (27) 

e- M '~ h) =l-2u (28) 

e° M = l-2u (29) 

where, 

u = (30) 



Consequently, 

[{2a-2b-\)u-a + b + l] 

Ci = 2 b y^2u) r r- (31) 



c 2 =- 



[(2a + 2b-l)u-a-b + l] 
2Zj(VT r 2^ 



a+6-1 



1 - 2m + (c/ + C 2 'r b ~ 2h )\ [(a + b)C l r b 2b +(a- b)C 2 J " +2 r 6 



(32) 
(33) 



For e"to be definitely positive in the region < r < r h , we must have 
Cj , C 2 ) .Thus in view of (3 1) and (32) we have; 



a+b-l a-b-l 

< u < 



2a + 2b-\ 2a-2b-\ 



(34) 



The central values of pressure and density are infinite, however, the limiting 

dp 



value of their ratio is finite and equal to the limiting value of 



dp 



f \ 
P 



^ dp^ 

\ d Pj r . >„ 



[2a + 2b-\ + b(n + 3>\a + b)] 
[l + b(n + 3\a + b)] 



(35) 



It has been calculated that for values of n in the interval [-1.3684, -1.058), the 
right hand side of (35) is negative thus making adiabatic sound speed imaginary and 
corresponding to the values of n in the interval [-1.058, -1], the right hand side of (35) 
is positive and less than 1 i.e. causality principle is obeyed at the centre ( Table- I). 
Hence for meaningful solutions n ranges within the interval [-1.058, -1]. 
Variation of central pressure- central density ratio / square of adiabatic sound speed at 
the centre and parameters of indices of e" for different values of n 



Table - 1 



S.N. 


n 


L 


a 


b 


P 1 \ dp \ 


1 


-1 


3.1667 


2.08335 


-0.58339 


0.333 


2 


-1.01 


3.1870 


2.07773 


-0.59714 


0.2656 


3 


-1.02 


3.2074 


2.07228 


-0.61058 


0.2037 


4 


-1.03 


3.228 


2.06699 


-0.62367 


0.1465 


5 


-1.04 


3.2486 


2.06184 


-0.63644 


0.0892 


6 


-1.05 


3.2693 


2.05680 


-0.64890 


0.04388 


7 


-1.058 


3.2860 


2.05294 


-0.65869 


0.00749 



5. Particular Members of Class 

In this section we shall present a detail study of the particular solutions corresponding 
to n = -1 and n = -1.05. 

(I) for n = - 1 , the solution is 

e"=( C ir °- 5 +C 2 r L667 f 

e A =- r — r + M(l.5C 1 +0.33C 2 r 117 ) 

[l.5C 1 r- 117 +2.667C 2 ] C, V ' 



~~4~ P 



Ar 



1.17 



0.38 



[l.5C ir - 117 +2.667C 2 ] + C 



(l.5Q 



+ 0.33C, r 



1.17 j 



(2Q +4.33C 2 r lA1 )_ ] 
(C, +C 2 r" 7 ) 



(-5.0004C! -5.7788C 2 r 117 )^ 
[l.SC^- 117 +2.667C 2 f 



272C r 1 ' 17 
U.Z/2C 2 r +Q43 

C, 



In view of (34) e v to be definitely positive in the region < r < r b , then 
0.25 < u < 0.3846 

Corresponding to u - 0.3 and in view of (31), (32) and (33) the constants are 
Qr, 05 =0.497, C 2 r h lM1 =0.1355, Ar b 2M =-0.2260 



In Table II the march of pressure, density, pressure - density ratio and square of 

adiabatic sound speed— is given for u = 0.3. We observe that pressure and density 

dp 

decrease monotonically with the increase of radial coordinate, pressure - density ratio 
and square of adiabatic sound speed is positive and less than 1 throughout within the 
ball. 



Table II 

The march of pressure, density, pressure - density ratio and square of adiabatic sound 
speed within the ball corresponding to n = -1 withw = 0.3 . 



X 


SnG 2 
4 Ph 

c 


%nG 2 
2 P r b 

c 




1 f dp\ 
c 2 [dp] 





oo 


oo 


0.333 


0.333 


0.1 


16.6018 


42.9419 


0.3866 


0.3708 


0.2 


4.7525 


10.9999 


0.4320 


0.4215 


0.3 


2.2907 


5.1544 


0.4444 


0.4801 


0.4 


1.3257 


3.1322 


0.4232 


0.5300 



0.5 


0.8273 


2.2045 


0.3752 


0.5960 


0.6 


0.5269 


1.7025 


0.3094 


0.6450 


0.7 


0.3204 


1.3990 


0.2290 


0.6933 


0.8 


0.1859 


1.1999 


0.1549 


0.8146 


0.9 


0.0807 


1.0609 


0.07560 


0.7905 


1.0 





0.9588 





0.7439 



(II) For n = -1 .05, the solution is 



Ar 



1.2978 



0.4197 



[l.4078C 1 r- 12978 +2.7057C 2 ]° 95 C, 



(l.3375C, +0.0397C 2 r 1 2978 ) 



SnG 



P = 



Ar 



1.2978 



0.4197 



[l.4078C 1 r- 12978 +2.7057C 2 f 5 C, 

(l.8175C, +4.411C 2 r 1 2978 ) , 

(c, +C 2 r L2978 ) 



(l.3375C, +0.0397C 2 



r 1.2978 j 



> X 



$7lG 



f(-4- 



97C, 



-6.217C 2 r 



1.2978 



11.95 



0.03 83C 2 r 



1.2978 



c, 



+ 0.439 



[l.4078C 1 r" 1 - 2978 +2.7057C 2 J 
In view of (34) e v to be definitely positive in the region < r < r h ,then 
0.2246<w< 0.3866 



Corresponding to u - 0.3 and in view of (31), (32) and (33) the constants are 



Qr, 04078 =0.4657, Qr/ ™ 57 =0.1667, Ar™ 2 =-0.2446 



In Table III the march of p,p, 



P 



and- 



dp 



within the fluid ball for u - 0.3 . 



pc~ c dp 

We observe that pressure and density decrease monotonically with increase of radial 



coordinate — ^- and — 
pc c 



1 dp^ 



and less than 1 throughout within the ball 



Table III 



r/ 
/ r b 


SnG 2 

4 Ph 

c 


%7tG 2 
2 P r b 

C 


P 

pc 2 


1 ( dp\ 

c 2 \dp) 





oo 


CO 


0.04388 


0.04388 


0.1 


4.9136 


44.2266 


0.1111 


0.0854 


0.2 


2.11629 


11.4429 


0.1849 


0.1370 


0.3 


1.2867 


5.3841 


0.2389 


0.2059 


0.4 


0.86436 


3.3462 


0.2583 


0.2842 


0.5 


0.59813 


2.3905 


0.2502 


0.3461 


0.6 


0.41072 


1.8728 


0.2193 


0.4375 


0.7 


0.27033 


1.5589 


0.1734 


0.5238 


0.8 


0.16126 


1.3459 


0.1198 


0.6428 


0.9 


0.07082 


1.2068 


0.0586 


0.6930 


1.0 





1.09915 





0.7350 



In Table IV, the variation of surface density (p b ), maximum Neutron star mass 
(multiple of solar mass M@), linear dimension corresponding to u - 0.3 with 
constant surface red shift Z b = [(1 - 2 5 -1] « 0.581 1 and different values of n . 

Table IV 



S.N. 


n 


%7lG 2 
2 Pbh 

C 


M/ 
/M & 


«2 r b Kms 


z h 


1 


-1 


0.9588 


3.24 


32.09 


0.5811 


2 


-1.01 


0.9918 


3.31 


32.64 


0.5811 


3 


-1.02 


1.029 


3.36 


33.26 


0.5811 


4 


-1.03 


1.0435 


3.39 


33.48 


0.5811 


5 


-1.04 


1.079 


3.44 


34.05 


0.5811 


6 


-1.05 


1.099 


3.48 


34.36 


0.5811 


7 


-1.058 


1.1372 


3.52 


34.69 


0.5811 



We observe that by decreasing the values of n the surface density, the mass and the 
linear dimension of Neutron star are increasing in nature. 



6. Conclusion 

For meaningful realistic fluid ball pressure and density must be positive and 



monotonically decreasing with the increase of radial coordinate r, and 4r 

pc c 



dp 
ydpj 

must be positive and less than 1 throughout within the ball. For meaningful solution 
e" = ( C x r a +C 2 r p f , where a and /? should be positive so that e v at centre 
becomes non-singular (positive and finite) but converse is not true. 
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